Interface thermal resistance between dissimilar anharmonic lattices 
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We study interface thermal resistance (ITR) in a system consisting of two dissimilar anharmonic 
lattices exemplified by Fermi-Pasta-Ulam (FPU) model and Frenkel-Kontorova (FK) model. It is 
found that the ITR is asymmetric, namely, it depends on how the temperature gradient is ap- 
plied. The dependence of the ITR on the coupling constant, temperature, temperature difference, 
and system size are studied. Possible applications in nanoscale heat management and control are 
discussed. 

PACS numbers: 67.40.Pm, 63.20.Ry, 66.70.+f, 44.10. +i, 
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When heat flows across an interface between two differ- 
ent materials (phases), there exists a temperature jump 
at the interface, from which we can define an interface 
thermal resistance (ITR): 
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where J is the heat flux density, namely, the heat flow 
across a unit area in unit time, AT is the temperature 
difference between two sides of the interface. This prob- 
lem has been caught attention as early as 1941 when 
KapitzaQ discovered the temperature jump at an inter- 
face between solid and liquid. Continuous efforts have 
been devoted to this problem (see Review article |2J, 
and the references therein) since then. More recently, 
the temperature jump between liquid and vapor has also 
been observed experimentally|j| and studied by com- 
puter simulations |4j. 

A general theory to describe the temperature jump at 
interface of different materials (in different phases) is still 
lacking. There are two popular theories for harmonic lat- 
tice: diffuse mismatch theory 2] and acoustic mismatch 
theory pj- The acoustic mismatch theory regards the two 
media as two elastic continua, while the diffuse mismatch 
theory assumes that at the interface all phonons are diffu- 
sively scattered. When anharmonic interaction is taken 
into account, superposition theorem fails, thus no ana- 
lytic theory can be worked out. However, as we shall 
show in the following that, if the anharmonic interac- 
tion is considered, more interesting phenomena arise, the 
problem becomes theoretically more challenging, and of 
course it is closer to reality because harmonic is just a 
first order approximation. On the other hand, as the 
rapid development of nanotechnology, low dimensional 
nano scale systems such as nanowires and nanotubes can 
be easily fabricated in the lab. At the nano scale, the sys- 
tems are of finite size, more precisely, they are discrete, 
and therefore, the continuous theory such as the acoustic 
mismatch theory will definitely not be suitable for such 
systems. The interface resistance in such nano systems 



becomes more and more important and has potential ap- 
plications in nanoscale heat control and management [(j. 

In this Letter, we study the ITR in the system com- 
prising two dissimilar anharmonic lattices. The system, 
illustrated in Fig. 1, consists of two chains of N oscilla- 
tors. The left part is a chain of N harmonic oscillators 
on a substrate whose interaction is represented by a si- 
nusoidal on-site potential. The right part is a chain of 
TV anharmonic oscillators. The two parts are connected 
by a spring of constant ki n t , and the Hamiltonian of the 
system is, 
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where Hfk is the Hamiltonian of the left part 
which is in fact the Frenkcl Kontorova (FK) model, 
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Hfpu is the Hamiltonian of the right part which is the 
Fermi-Pasta-Ulam (FPU) model, Hfpu — Si=jv+i 2^1 + 
kFPU \{xi+i - Xi - a) 2 + | (x i+ i - Xi - a) 4 . 
Fixed boundaries are used, i.e. xo = and 
x 2N +i = Nb + {N + l)a. The FPU model is a 
representative anharmonic lattice without on-site poten- 
tial and the FK model is the one with on-site potential. 
Both models have been widely used to study different 
problems in condensed matter physics and nonlinear 
dynamics0,0j. In particular, the FPU model has played 
an important role in the development of computational 
physics and nonlinear dynamics 0- Recent years have 
witnessed an increasing interests in the study of heat 
conduction with these two models |9i llfj. llll Il2l 1 1 3| . 

In our simulation, Nose-Hoover heat baths are used. 
The system parameters are kpK = L, j3 = 1, a = b = 1, 
m = 1, and V = 5. Tj, = T (l + A),T R = T (l - A). 
kFPU = 0.2 unless otherwise stated. The local tempera- 
ture at site n is defined as T n = m{x n 2 ), where ( ) stands 
for temporal average. We calculate the physical quanti- 
ties after a time that is long enough to allow the system 
to reach a non-equilibrium steady state where the local 
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FIG. 1: Configuration of the system. The left part is a 
Frenkel-Kontorova model and the right one is a Fermi-Pasta- 
Ulam model. The two parts are connected by a spring of 
constant ki nt . The two particles on the left and right ends 
are put into contact with heat baths of temperature Tl and 
Ta, respectively. 
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FIG. 2: (a) T/T versus lattice site for T = 0.07, 0.09, and 
0.3. |A| = 0.5. (b) (T - To)/(T |A|) versus lattice site, for 
different |A| = 0.2,0.4, and 0.7 with fixed To = 0.09. The 
solid symbols are for the cases of A > 0, and the empty ones 
are for the cases of A < 0. In both (a) and (b) N = 50. 



heat flux is constant along the chain. 

Figl^a) shows the normalized temperature, T/Tq, 
along the lattice site for three different temperatures, 
T o =0.07, 0.09, and 0.3, with fixed |A| = 0.5. Fig. □» 
shows (T - r )/(T |A|) versus lattice site for |A|=0.2, 
0.4, and 0.7 with fixed To = 0.09. It is obvious that, in 
all cases, there exists a temperature jump (discontinuity) 
at the interface. The jump depends not only on tempera- 
ture T but also the temperature difference A. The most 
interesting thing is that, the jump is asymmetric, namely, 
it depends on whether A is positive or negative. 



10° 
10 5 

io 4 
io 5 
io 2 

0.01 



(a) cx 



o x R 
oo - 

o 



\ 
o 



R 



o 
o 

6 10 

Q 



▲ A 

/ £ 

A - A -J 

A^ 



o 10 

,0.01 

CX 



0.1 T 



•» s sag 



0.1 



io 3 



io 4 



10' 



io z 



(b) ,o' 

io 5 



10 4 

io- 1 



10 

ur 
io ! 
ltr 



O 



R 



o 



R / 



O 



o 



|AT| 0.1 



-O' 
-•- 



,o 



R 



o.i 



|A| 



1 



FIG. 3: (a) R± versus To. |A| = 0.5. Inset is heat current J± 
versus T . (b) R ± versus |A| for T = 0.12. Inset of (b) is R± 
versus interface temperature Tint = (T„ 4 + T^ t )/2, and R± 
versus interface temperature jump, AT. In all cases iV = 50. 



The asymmetric feature is well captured by ITR shown 
in Fig. □ R± denotes the ITR when A > and A < 0, 
respectively. Fi^3Ja) shows that, at intermediate temper- 
ature, To <= (0.02,0.1), the ratio, R-/R+, changes from 
about 10 to 1000. Both R + and R- decrease with To un- 
til To « 0.2 and then become approximately constants. 
The dependence of R± can be explained from the tem- 
perature dependence of the phonon spectra, as we shall 
discuss fully later on. In the inset of FigJSfa), we show 
the heat currents J± versus temperature To. J± is the 
heat flux for A > and A < 0, respectively. Heat cur- 
rents increase very fast (in several orders of magnitude) 
as To is increased over two orders of magnitude. 

In Fig. Eft)), we draw R± versus |A| for T = 0.12. 
It shows that as |A| increases, R+ does not change too 
much, it is always in the order of 200, whereas i?_ in- 
creases more than three orders of magnitude. In the 
case of |A| = 0.8, R-/R+ can be larger than 1,000 and 
\J+\ w 2,00Q|J_|. The dependence of ITR on the in- 
terface temperature, T„ t and the interface temperature 



jump, AT 
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Tj„ t are shown in the inset of Fig. 



EJb). Where T in ' t is the temperature for the particle on 



3 



the left /right side of the interface, respectively. 

An obvious conclusion from above results is that the 
ITR between dissimilar anharmonic lattices is asymmet- 
ric. 

To understand the physical mechanism of this phe- 
nomenon, we need to invoke the energy band theory. 
However, due to the presence of anharmonicity, an an- 
alytic approach seems impossible. We shall rather take a 
qualitative approach, which can also give us useful infor- 
mation in two extreme cases, namely, low temperature 
limit (regime) and high temperature limit (regime). 

FK model. At very low temperature, the particles 
are confined in the valleys of the on-site potential. By 
linearizing equations of motion one can easily obtain the 
phonon band, yfV < oj fk < y/V + Akpx- On the other 
hand, in the high temperature limit, the particles have 
large enough kinetic energies to jump out the valleys. 
The on-site potential becomes negligible, the FK model 
degenerates to an harmonic one, < 0Jp K < 2^/kpK- 

FPU model. There exists a threshold temperature 
Tthfa 0.1)|l0j, below which the FPU model becomes a 
harmonic one, one thus has < Up PU < 2^/kppjj. When 
To >> Tth, the anharmonic term is dominant. In this 
regime, a rough theoretical estimate yields < ojp PU < 
C (Tk F pu(3) 1/4 with Co = 2V2^r(3/4)3 1 / 4 /r(l/4) « 
2.23, where T is the Gamma function. 

The spectra of the interface particles are shown in 
Fig'EI a ) an£ l 3 ) f° r A > and A < 0, respectively, 
and compared with the above analytical analysis (the 
shadow regions). In the first case, when Tl = 0.15 and 
Tr = 0.01, we can approximately regard the FK lattice 
as at high temperature limit, and the FPU lattice at low 
temperature regime. With kpK — 1 and kppjj — 0.2, we 
have, u>p K € [0, 2], and Up PU € [0, 0.89], which are quite 
close to the numerical ones. 

On the other hand, when the temperatures of the 
two thermal baths are swapped. Approximately, the 
FK lattice is at low temperature limit, and the FPU 
is at high temperature regime. In this case, according 
to above analysis, the phonon band of the FK model 
is u>^ K £ [2.24,3.00] for V = 5 and k FK = 1, and 
uptpu € [0,0.94] with T rs 0.15 and (3 = 1 for the FPU 
model. They are also very close to the numerical results 
shown in Fig^Jb). 

It is clear from FigQ] that, in the case of Tl > Tr, 
the phonon band of the FK lattice overlaps that one of 
the FPU lattice, therefore, heat can easily flow from the 
FK lattice to the FPU lattice as is demonstrated in Fig|3| 
by small R+ and large J+. Conversely, when Tl < Tp, 
by appropriately chosen parameters, a large gap between 
the phonon bands of the FK lattice and the FPU lattice 
can be formed (see FigQJb)), which inhibits heat flow 
from the FPU lattice to the FK lattice as is manifested 
by a large i?_ and small J_ in Fig|2| However, as To 
increases, the gap becomes narrower and narrower, and 
eventually disappears when Tq surpasses a certain value. 
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FIG. 4: Phonon spectra of the two particles at interface and 
schematic phonon bands for the FK and the FPU models in 
two extremes, (a) T L = 0.15 and T R = 0.01. (b) T L = 
0.01 and Tr — 0.15. The shadow regions are the analytical 
estimates for the FK and FPU models, (see text for more 
explanation.) 

Indeed, when To is large enough the two phonon bands 
overlap leading to a constant i?_ as is seen in Fig 3^a). 

In order to quantify above asymmetric ITR, heat cur- 
rents, and find the relationships with the overlap of the 
phonon bands of the two lattices, we introduce, 

S± = Jo^Pl^W^ (3) 
Jo" P l ( UJ ) duJ Jo" P r{u)dLO ' 

S± correspondes to the overlap for A > and A < 0, re- 
spectively. Note that J °° Pi >r duj = T^f . In Fig|3]we plot 
R—/R+ versus S+/S- and in the inset we show |J + /J_| 
versus S+/S-. Best fit gives R-/R+ - {S+/S-) Sr 
with 5 R = 1.68 ± 0.08, and \J+/J-\ ~ {S+/S-) Sj , with 
5 j = 1.62 ±0.10. This picture indeed illustrates that the 
ITR and heat current correlate strongly with the overlap 
of the phonon spectra of the two interface particles. 

As the system contains many adjustable parameters, it 
is worth investigating the dependence of the ITR on these 
parameters. Fig. HJa) is R± versus hint, which shows 
both R± decreases with increasing ki n t ■ This means that 
strong coupling favorites heat transport. Fig|SJb) shows 
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FIG. 5: R-/R+ versus S+/S-. Inset, J+/J- j versus S+/S- 
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FIG. 6: (a) R± versus the coupling constant Aw. (b) R± 
versus k,FPu- (c)R± versus the system size, N. For all three 
cases, T = 0.07 and |A| = 0.5. N = 50 for (a) and (b). 
kint = 0.05 for (b) and (c). 



R± versus kppu and FigEIc) demonstrates the finite 
size effect. R + increases slightly whereas R- decreases 
slightly with TV. This can be understood from the fact 
that with fixed T + and T_, when A > 0, the larger N 
the smaller the T^ nt , thus the larger R + ; conversely, when 
A < 0, the larger N, the larger T^ nt thus the smaller 
The asymmetry of the ITR and the heat current in the 



model studied is due to the anharmonicity (nonlinearity) 
of the two lattices. In the low temperature limit, such 
as T < 0.01, in which both the FPU and the FK mod- 
els can be approximated by harmonic lattices, the asym- 
metry property will vanish because the phonon band of 
harmonic lattice is temperature independent. On the 
other hand, in the high temperature regime, in particular, 
when T >> V/(2n) 2 such that the on-site potential can 
be neglected, then the phonon band of the FK model be- 
comes temperature independent. In this case, the asym- 
metry effect becomes minimal as we can see from Fig. 
3 when To > 0.2. However, the asymmetry will never 
disappear as long as the anharmonic term in the FPU 
model is still present. 

In summary, we have studied the ITR between two dis- 
similar anharmonic lattices. It is found that the ITR is 
asymmetric, which is believed to be very general as the 
two anharmonic lattices used in this Letter are two repre- 
sentative ones widely studied in different fields of physics. 
The asymmetric property might be useful in heat control 
and management 6J. In particular, the very high ITR 
might find applications in building thermal insulator. 
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